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Abstract 



From an estimate of the meson-meson inelastic scattering at 5 GeV it is 
concluded that a typical strong phase in B decays to two mesons is of order of 
20°. For a particular final state an estimate of the phase depends on whether 
that state is more or less probable as a final state compared to those states 
to which it is connected by the strong interaction S matrix. 
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I. INTRODUCTION 



Strong final state interactions play an important role in the analysis of CP-violating 
effects in B decays. Direct CP violation such as a difference in rates for B + — > F and 
B~ — ► F vanishes in the limit that there are no strong phase shifts. Final state phases play 
a critical role in amplitude analyses of a set of B° decay experimental results. 

An approach to final state phases in inclusive decays was given by Bander, Silverman, 
and Soni [|IJ. For decays corresponding to the transition b — > uus they considered at the 
quark level b — *> ccs — > uus where the second transition on mass shell yielded the final state 
phase. Whether or not this is reasonable for inclusive decays, its application to exclusive 
decays such as B — > Kir has been criticized [Q because the major final state interactions of 
Kit are "soft" scattering to Kutc and not to cc states. 

We concentrate here on the decay of B to two mesons, referring to B — > tttt to be specific. 
Arguments have been given that the final state phase shifts should be small. For example, 
Bjorken [|J argues that there is little final state scattering in B — > tttt because the B decays 
directly to colorless qq pairs that do not interact as they evolve into tttt. Taken literally 
this is not correct since the s-state scattering at 5 GeV is expected to be sizable |5J as we 
discuss below. In the present note we seek to analyze the relations between the weak decay 
amplitude and the strong interaction S-matrix of final states that might lead to large or 
small final state phases. 



II. MULTICHANNEL FINAL STATE INTERACTION 

Consider the decay matrix element (f out \Oi\B) for the B meson into the hadronic final 
state /, where Oi is a weak decay operator. The strong interaction S-matrix is defined with 
the "in" and "out" states by 

s ff > = (rvn- (i) 

We choose states that are eigenstates of J not of individual meson momenta. The phases of 
the in and out states are fixed by the time reversal transformation T: 

T\r in ) - <m (2) 

With this phase convention, time reversal invariance of strong interactions requires that Sff> 
be symmetric matrix: 

s ff = s f'f- ( 3 ) 
Applying time-reversal operation on Mf = (f out \Oi\B) , one obtains 

M f ^ (B\TOiT- l \f in ). (4) 

If one inserts a complete set of out states and uses Eq.(^|), this relation becomes Mf = 
J2f> SffiMfi for a T-even decay operator Oi. One can express it in the operator form as 
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M = SM*, 



(5) 



where M is represented in a column vector. This matrix equation is formally solvable as 

M = S 1/2 M°, (6) 

where M° is an arbitrary real vector of the same dimension as M. If one uses the eigenstates 
| a) of the S-matrix as a basis, Eq.(|6|) reduces to the Watson theorem: M a = M®e lSa . We 
thus may consider the vector M° as representing the decay amplitude in the absence of the 
final state phases due to the strong interaction. Since M and M are related by a unitary 
matrix, it holds that £/ \M f \ 2 = £/ |Mj?| 2 . 

If one subtracts the complex conjugate of M from both sides in Eq.(|) and divides by 2i, 
the familiar form ImM = tM* emerges for the imaginary part of M, where t = (S — l)/2i. 
In components, it reads 

lmM f = J2t fr M},. (7) 
/' 

This form is commonly derived starting with Lehmann-Symanzik-Zimmermann's reduction 
formula. In applications of interest, the weak decay Hamiltonian H w is given in the form 

i 

where Aj is a combination of the CKM matrix elements and Oi is a T-even operator. It is 
to be understood that Mf is to be evaluated separately for different operators (9;. 



III. STRONG INTERACTION S-MATRIX 

When two mesons such as tt + tc~ interact in the s-state, we believe that they will scatter 
into a large number of multi-particle final states. Indeed we expect similar inelastic behavior 
for all partial waves of I < kr where r is a characteristic hadron radius. The sum over 
these partial waves can be described by a diffractive scattering formula such as that given 
by Pomeron exchange. For the case of meson-meson scattering we extrapolate from the 
analysis of meson-baryon and baryon-baryon scattering and write for the invariant elastic 
scattering amplitude , 

T(s,t) =ia tot se u , (9) 

where the constant in front has been fixed by the optical theorem. Defining the Pomeron 
residue by (3{t) = <Jtote bt , we obtain with the factorization relation (3(t)MM'P(t)pp = 

(3(t)MpP(t)M'p Q 

a™~12mb, < o f~10mb, ff**~8mb, (10) 



The factorization can be proved only for simple Z-plane singularities. It is an assumption for the 
Pomeron. 
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where a P ot = 37mb, a^t — 21mb, and af Q f = 17mb || have been used for the diffractive 
contribution of a tot at y/s — m B- For the diffractive peak width, the factorization gives 

b™ ~ 3.6GeV~ 2 , b K7T ~ 2.8GeV" 2 , ~ 2.0GeV~ 2 (11) 

if we choose b pp ~ 5.0GeV" 2 , ~ 4.3GeV~ 2 , and b Kp ~ 3.5GeV~ 2 0. For Dvr scattering, 
we use the quark counting rule for a to t and the assumption that the charmed quark interacts 
with the light quarks much more weakly. Then we obtain a crude estimate 

°tot ~ 2 to *' ^ ^ 

and 6 D?r is a little smaller than 6^^. 

Projecting out the s-wave from the amplitude in Eq.(|9]), 

a z=0 (s) = /° T{s,t)dt (13) 

1D7TS J— s 

yields 

{0.16 (tttt) 
0.17 (Kir) , 
0.18 (KK) (14) 
0.12 (£)7r) 

at yfe = 5 ~ 6GeV. 

Extraction of the s-wave amplitude from the diffractive formula may arouse suspicion 
since one thinks of diffraction as a peripheral process. || It would be better to consider 
T(s, t) as describing the scattering from an absorbing gray sphere of radius r. The values 
of ai up to / ~ kr vary gradually with I thus adding up to a large forward peak. As a 
result about 90% of the contribution to the integral in Eq . (|13|) comes from \t\ < lGeV 2 even 
though the I = amplitude itself is, of course, independent of t. 

In what follows we use the estimate 

Sff * 0.7 (15) 

corresponding to 

a l=0 = t ff = = 0.15z. (16) 

This corresponds to the case of a gray sphere with an inelasticity of 0.85. In the limiting 
case of a black sphere Sff goes to zero and the inelasticity goes to 0.5. If one goes beyond 
the diffractive scattering approximation, Sff is not purely real. In the Regge theory, the 
real part arises from exchange of the non-Pomeron trajectories such as p and fa. In n^p 
scattering the real-to-imaginary ratio of 10 ~ 20% was observed in the forward scattering 
amplitudes at ^/s = 5 ~ 6 GeV. [0] We can make an estimate of the real part for meson- 
meson scattering by using the factorization. We first determine the Regge parameters at 
t — from the total cross section differences || and then extract their t-dependences from the 
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angular dependence of the differential cross sections. [0J The analysis is simpler if exchange 
degeneracy is imposed. The smaller a to t and the larger p-f 2 residues tend to enhance the 
real to imaginary amplitude ratio for inr scattering over irp scattering, while the smaller 
b 71 ^ partially compensates the trend. Particularly for 7r + 7r°, the real parts of the p and f 2 
terms add up close to 30% of the imaginary part. However, our major goal is to understand 
the implications of the sizable inelastic scattering; for this purpose we use the simplifying 
approximation that Sff is real. 



IV. TWO CHANNELS 

The relation (|5|) was studied in the case of two channels U assuming that the diagonal 
S-matrix elements Sff are purely real. This requirement on the S-matrix turns out to be so 
strong in the case of two channels that there is only a single parameter left: 



S=[ f f* Tn Ae) ^ f , (17) 




/ cos 20 isin20\ 
" = \ isin26 cos 26 ) ' 

When S 1 / 2 is computed from S and substituted in Eq.(||), a simple relation results: 

M x = M i cos 9 + iM 02 sin 9 , (18) 
M 2 = iM 01 sin 6 + M 02 cos 6. 

The phase of the decay amplitude in the channel 1 is large if the particle decays preferentially 
to the channel 2, while it is small if the channel 1 is the dominant decay channel for small 
values of 9. Though it is an interesting conclusion, this picture turns out to be specific to 
the two-channel case. When we add one more channel in the final state, there are three 
real parameters even after \mSff = is assumed. The nice simple relation of Eq.(18) does 
not hold any longer. If we go to iV final channels, there are N(N — l)/2 real parameters 
even with ImS// = 0, and no meaningful prediction results. Therefore we must change our 
strategy in studying the case of N 3> 1 such as the B decay. 



V. RANDOMNESS OF WEAK DECAY AMPLITUDES 

In the presence of many decay channels, strong interactions are so complicated that it 
is beyond our ability to predict final state interactions accurately. We must substitute lack 
of our knowledge with reasonable dynamical assumptions and/or approximations. In search 
of such an assumption, we notice that since Mf and Sff come from two different sources, 
weak and strong interactions, the phase of product SffM% for f ^ f takes equally likely a 
positive or a negative value as /' is varied with / fixed. While Mf/ is related to Mf (f^f) 
by rescattering, there exist so many states that the influence of / on /' can be disregarded. 
We therefore introduce the postulate that the phase of SffM% takes random values as /' 
varies. It should be noted that randomness is postulated here for the relative phase and 
sign of the decay matrix element to the S-matrix element, not for the dynamical phases and 
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mixing of S-matrix as it was introduced in the random S-matrix theory of nuclear physics. 
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We start our analysis with Eqs.(|7]) and (fKI). Taking the /' = / term in the sum to the 
left-hand side in Eq.(^) and using £// — ilmtff, we write Eq.(^) in the form 

(1 + it ff )lmM f - t ff ReM f = £ t ff M* fl . (19) 

f'+f 

The first and the second term of the left-hand side are real and imaginary, respectively, for 
Ret// = 0. Given the estimate of Eq.flTo]) the coefficient of the first term is much larger than 
that of the second term and we consider this primarily as an equation for ImMj. With the 
randomness postulate, the phase of Mf is equally often positive or negative if we consider 
a large ensemble of final states /. It is some kind of average of the magnitude of the phase 
of Mf, not values for individual Mf, that we can study with our randomness postulate. 
For this purpose, we take the absolute value squared for both sides of Eq. (|I9"D . Then the 
right-hand side is: 

RHS = jJ2 S ff' M *f> S *ff" M f", (2°) 

where tfy = Sff/2i has been used for /' ^ f. The random phase postulate allows us to 
retain only the terms of /' = /" and to reduce the double sum to a single sum: 

RES^Ws f fS} lf \M r \ (21) 



- A \M}\ £ 1%'f, (22) 



where the second line defines |M?,| as the weighted average of the decay amplitudes into 
states /'. Then, using the unitarity of S-matrix, we reach 

1 



RHS~ -{1 - Sj f )\M},\. (23) 



While our estimate of Sff is made on the basis of the Pomeron exchange, it should be noted 
that contributions to Sff from quantum number exchange may be important in determining 
\Mf,\ from Eq. (|22l) if they correspond to states /' with large values of |M|, |. Identifying 
Eq.(p3"D with the absolute value squared of the left-hand side of Eq. (p!9[) , we obtain the 
prediction of our random phase approximation: 



(1 + Sf f f{lmMff + (1 - 5 // ) 2 (ReM / ) 2 = (1 - Sj f )\Mf,\. (24) 

Defining p by 

p = ]Mj\ 1/2 /\M f \, (25) 
\M f \ 2 = (ImM/) 2 + (ReM/) 2 , 

the ratio of the imaginary-to-real part of Mf is solved from Eq. fl24]) as 
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(ReM f y 



= tan 2 Sf 



r 2 ( P 2 



T 



p 2 T 2 



where 



T 



l + S 



ff 



(26) 



(27) 



Note that r 2 is equal to the ratio of elastic to inelastic scattering cross section a e i/a ine i of 
the relevant partial wave. Since the left-hand side of Eq.(^) is nonnegative, r and p are 
constrained for Sff > by 



r 2 <p 2 < 1/r 2 



For S 



ff 



0.7, 



T 



0.18, 



\M f/ \ 2 and 



(28) 



(29) 

M/| 2 to differ too 



so that rescattering among the final states does not allow 
greatly in magnitude. In the weak limit of rescattering (r — ► 0), of course, Eq.(^) allows 
any value for p. In the black sphere limit (r — ► 1) Eq.(|26|) is useless and Eq.(^) constrains 
p = 1. Our approach is only useful to the extent that inelastic scattering dominates very 
much over elastic for the final state /. It should be noticed that Eq.(|26|) reduces to the 
two-channel case Eq.(18) with 



M 2 ° 



Mi 



P 



T 



p2 T 2 



(30) 



Our random approximation amounts to lumping all inelastic channels together as if they were 
a single inelastic channel with an "average" decay amplitude. However, we now interpret 
this as something like the standard deviation of the phase for an ensemble of independent 
final states / with a given value of p. 

If the relevant states /' were similar to the state /, then we might expect p to be close 
to unity. For p = 1 Eq.(p6|) reduces to 



tan 2 5 



1-S 



f 



T 



ff 



l + S 



ff 



(31) 



sin 8j 



1-S 



ff 



With Sff = 0.7 this gives \5f \ ~ 23°. Thus a typical value of the final strong interaction phase 
in this case is not small. This result for a typical state has a simple heuristic interpretation. 
The original real decay amplitude M° is reduced as a result of absorption by a factor a, but 
an imaginary term arises due to rescattering from other states. Since the total decay rate is 
not changed by final-state scattering the final value of \Mf\ for a typical state will be equal 
to \M?\. Thus M f takes the form 



M f = M°[a + iVl-a 2 



(32) 



ImMf 
ReM f 



V 1 — a 2 1 a. 
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This agrees with the result above if the absorption factor is identified as 

a = + S ff )/2. (33) 

Any argument that a final state phase is small must be an argument that p is small. It 
should be noted that p depends on the particular final state / and on the weak interaction 
operator O^. The quantity \Mj,\ is an average of the square of the decay amplitude to state 
/' via Oi weighted by the square of the scattering amplitude from / to /' (cf Eqs.(p2|) and 
(Iffi)). Thus a value of p much smaller than unity means that on average the states to which 
/ scatters are much less likely than / to be final states in the decay due to operator Oi. 
Conversely, if / is a particularly unfavored final state p may well be above unity. Figure 1 
shows the dependence of the phase on p for the choice of Sff — 0.7. 

VI. POSSIBILITY OF SMALL OR LARGE STRONG PHASES FOR TWO-BODY 

CHANNELS 

While our randomness postulate leads to sizable strong phases for "typical" decay chan- 
nels, dynamical arguments are necessary to estimate the parameter p for a specific channel 
and a specific decay operator. We ask what dynamical arguments could let our prediction 
agree with Bjorken's argument which favors small phases for decays such as B° — > tt + tt~. 
The observed branching fraction of B° — > 7r + 7r~ is about 10~ 3 of the inclusive branching 
fraction for b — > dqq within large uncertainties. This might seem to indicate that 7r + 7r~ 
is not a favored final state and that the strong phase of 7t + 7t~ might be of order of 20°. 
However this conclusion can be evaded if certain conditions are satisfied. 

It should be noted that we are only interested in whether ti + it~ is a favored decay channel 
relative to those to which it is connected by Sff. Most of the states /' are multi-meson 
states with little jet-like character. It could be argued that these states are not favored as 
B decay products because they are not likely to develop from three quark jets into which B 
naturally decays. It is not clear whether this distinction is really operative for the energy 
m B . 

If the above argument were true, it could be considered as an interpretation of the 
Bjorken argument. In order to produce a tttt final state the final quarks must emerge as 
colorless pairs; alternative quark configurations are unlikely to hadronize into tttt. Thus p 
would be close to its minimum value and ImMj would be small. It is not true, however, 
that the final tttt state has little interaction, but the effect would primarily be a moderate 
absorption correction to the real part. In the two-channel reduction this corresponds to 
M 2 o/M 10 close to zero and the change of the real part of Mi from M 10 to Mxqcos^ would 
be considered as the absorption correction. 

An alternative possibility in which one might expect a large final state phase shift has 
been emphasized in a number of recent papers (T^|. These are cases in which Mf vanishes 
in the naive factorization approximation. An example is the "tree" amplitude proportional 
to X u = V u bV* s for the decay B~ — > K n~ . In this case the major contribution to Mf is 
expected to arise from rescattering from favored states so that p might be larger than unity. 
It can be argued |HJ that in this case, even though the final state scattering mainly goes to 
multi-particle states, the main contribution to the strong phase arises from the scattering 
of quantum number exchange involving two particle to two particle transitions 
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VII. CONCLUSION 



For a "typical" final state the strong final-state phase shift is not small; a typical mag- 
nitude is 20°. We understand this to be the magnitude averaged over the states that are 
interconnected by the final-state S'-matrix. We expect there to be sizable fluctuations about 
the average; in fact, since we cannot predict the sign of the phases, our analysis suggests 
that the algebraic average phase may be zero. 

A simple heuristic understanding of this phase is that the final state absorption reduces 
the value of the original real decay amplitudes whereas rescattering from other states pro- 
vides an imaginary amplitude. For a "typical" state the absolute value of the amplitude 
is not changed since the final state interaction does not change the total decay rate. Our 
magnitude estimate of ~ 20° is derived from the expected inelasticity of the meson-meson 
scattering. 

It may be possible to argue for a particular final state / that the phase is small. Any such 
argument must show that the states to which / is connected by the S'-matrix are generally 
less likely to be a decay product of B than /. Thus it might be argued that the four-quark 
operator leads easily to 7r + 7r~ whereas there are many states to which 7r + 7r~ scatters that 
are not easily reached directly via B decay. Conversely it might be argued that for states 
which are not easily reached via the four-quark operator, the strong phase is large. 
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FIGURES 



delta (degrees) 




FIG. 1. The strong phase 8f defined in Eq.(26) plotted against the ratio p for Sjf = 0.7. 5f is 
chosen between 0° and 90°. 
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